Scattering of periodic solitons 
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Abstract 

With the help of numerical simulations we study A^-soliton scattering (A^=3,4) 
in the (2+l)-dimensional CP^ model with periodic boundary conditions. When 
OO ■ the solitons are scattered from symmetrical configurations the scattering angles 

observed agree with the earlier n/N predictions based on the model on 3^2 with 
I standard boundary conditions. When the initial configurations are not symmet- 

ric the angles are different from n/N. We present an explanation of our observed 
patterns based on a properly formulated geodesic approximation. 
Oh! 

'f: : 1 Introduction 

> . 

^ , Physics in (2+1) dimensions is an area of much active research, covering topics 

?H I that include Heisenberg ferromagnets, the quantum Hall effect, superconductiv- 

ity, nematic crystals, topological fluids, vortices and solitary waves ilj. Most of 
these systems are non-linear. In their mathematical description the well-known 
family of sigma models plays a starring role. The simplest Lorentz-covariant soli- 
ton model in (2+1) dimensions is the CP^ or non-linear 0(3) sigma model. Its 
solutions, sometimes called 'lumps', are realisations of harmonic maps, a long- 
established area of research in pure mathematics. However, analytical 0(3) soli- 
tons solutions have only been found for the static case; the full time-dependent 
model must be studied using numerical methods and/or other approximation 
procedures 
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Sigma models are also useful as low dimensional analogues of field theories in 
higher dimensions. In effect, the 0(3) model in two dimensions exhibits confor- 
mal invariance, spontaneous symmetry breaking, asymptotic freedom and topo- 
logical solitons, properties similar to those present in a number of important field 
theories in (3+1) dimensions -like the Skyrme model of nuclear physics 3 . 

We are concerned with the planar CP^ model (both in its original and 
Skyrme-like versions) with periodic boundary conditions where the solitons are 
harmonic maps T2 82- A rich diversity of phenomena has been found in this 
model [H ISl El ) going beyond the two-soliton and annular structures one might 
expect by analogy with the model with standard boundary conditions, where the 
soliton fields are harmonic maps 9?2 U {00} « 52 1— > 

For planar systems, identical lumps initially placed at the vertices of a 
rectangular A^-gon, fired with equal speed to collide head-on at the centre of 
the polygon, scatter and emerge on the vertices of the dual polygon. Such tt/N 
scattering, studied for the CP^ model on 9^2 in reference 0, may be understood 
on symmetry grounds: the initial data has (dihedral group) symmetry and 
the time evolution respects that. But imposing periodic boundary conditions 
breaks the foresaid symmetry, and the interesting question whether dual-polygon 
scattering still holds for this case must be investigated. 

Through numerical simulations on T2 we confirm the tt/N scattering for 3 
and 4 identical lumps. The case N = 2 has being considered elsewhere |1] and 
it conforms to the well-documented scattering at 90°. We also look into non- 
symmetrical configurations (they do not scatter at tt/N) and explain the results 
using the geodesic approximation. 

In the following section we introduce our periodic CP^ model. The numerical 
procedure is explained in section |3J Section |1] analyses collisions between three 
solitons, and the case A^ = 4 is considered in sectional In section IHl we present 
our version of the geodesic approximation which includes the case of initial con- 
figurations that are not symmetrical. We find that our predictions (based on this 
approximation) are in agreement with what is observed in numerical simulations. 
We close with some concluding remarks in section [3 



2 The CP^ model on the torus 

The non-linear 0(3) model involves three real scalar fields (p{x^) = {4>i, 4'2, 4'3) 
which satisfy the constraint that V = {x^,x^,x'^) = {t,x,y) G T2 the fields lie 
on the unit sphere ^2: 

$.^=1. (1) 
Subject to this constrain, the Lagrangian density of the system is 

C = \id,$).{d'^$). (2) 

The invariance of the model under global 0(3) rotations in field-space is apparent. 
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The model is conveniently recast in terms of one independent complex field 
W (the CP^ formulation) related to (p via 

j:_r W + w . -w + w \w\ 



The Lagrangian ^ now reads 




. ^ \dtW\^-2\d,W\^-2\d,W\ 
(1 + |H^|2)2 ' 



where z = x + iy and z = x — iy. For all t, the fields W are mappings from the 
torus T2 to the sphere 5*2, ie, they satisfy the periodic boundary conditions 

W{z + mL + inL) =W{z), (5) 

where m,n = 0, 1, 2, . . . and the period L denotes the size of the square torus. 

The static soliton solutions of the model are thus doubly periodic functions of 
z, that is, elliptic functions that may be expressed through Weierstrass cj-function 
as [HI El 

^=^n^p|4> E%=Ev (6) 

J=l^i^ ^3) j=l j=l 

The complex number A is related to the overall size of the solitons, the zeros aj 
and poles hj determine their sizes and positions on T2, and the positive integer 
K is the order of the elliptic function W . 

The static energy density (or potential energy density) associated with © 
can be read-off from 

^ = ^ {i + \w\^Y ■ 

Pictures of this energy distribution reveal the familiar CP^ lumps, as those of 
figure n for example. 

The energy is related to the topological charge by the Bogomolnyi bound 

E > 2Tr\N\. (8) 

The instanton solutions correspond to the equality in solutions carrying 
N > {N < 0) imply dzW = {dzW = 0), the Cauchy-Riemann conditions for 
W being an analytic function of z (z). Note that the simplest non-trivial elliptic 
function is of order two, hence there are no single-soliton solutions on T2. 

We utilise the pure model (jlj) for = 3, but for = 4 the energy involved 
is so large that the well-known instability of the planar model (recall its confor- 
mal invariance means that the lumps can have any width) breaks the numerical 
procedure fast. A stabilising Skyrme term must be introduced for this case. 



3 



3 Numerical procedure 



In this paper we want to discuss time dependent solutions of the model ((3} 
and its Skyrme version [see equation where the time dependence describes 

the movement of solitons. As we do not have analytical time-dependent solu- 
tions, we resort to numerical simulations. For this we take fields of the form 
dnj as the initial conditions. Since during the simulations the field W may 
become arbitrarily large, we have preferred to run our simulations in the (p- 
formulation of the model. One returns to the CP^ formulation inverting @: 

W={l-(j)3)/{cj)i+i^2)- 

Strictly speaking, truly independent solitons can only be obtained in the 
asymptotic regime of large soliton separation, which really never happens on a 
compact manifold such as T2. However, each factor Wj = a{z — aj)/a{z — bj) in 
© (when aj ~ bj) roughly represents one soliton, providing a setting to studying 
more or less independent structures. The present work is limited to systems in 
the topological classes 3 and 4. These systems move, collide and scatter off upon 
being set into motion by boosting each Wj separately. The initial-value problem 
is then completely specified by giving both W{t) and dtW{t) at the initial time 
t = 0. 

For a square torus we have the so-called lemniscatic case |lUj where a posseses 
the simple Laurent expansion 
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Gj = Gj{L) G K. 



(9) 



For the accuracy of our calculations it is sufficient to compute the series 
to G5 (our coefficients for j > 5 are negligible): 



up 



Co 


= 1 


Gi 


= -0.7878030 


G2 


= -0.221654845 


G3 


= 9.36193 X 10^3 


G4 


= 7.20830 X 10"^ 


G5 


= 2.37710 X 10"^ 



We have employed the fourth-order Runge-Kutta method and approximated the 
spatial derivatives by finite differences. The Laplacian has been evaluated using 
the standard nine-point formula and, to further check our results, a 13-point 
recipe has also been used. Respectively, the Laplacians are: 



1 4 1 
4 -20 4 
1 4 1 

6 X {6xy 



1 

-1 12 



-1 

12 1 

-48 12 
1 12 1 

-1 
10 X ((5x)2 



(10) 
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The discrete model has been evolved on a Ux x Uy = 200 x 200 square periodic 
lattice with spatial and time steps 6x=5y=0.02 and 5t=0.005, respectively. The 
vertices of the fundamental period cell we have used for our simulations were at 

(0,0), (0,L), (L,L), (L,0), L = nxx6x = 4. (11) 

Unavoidable round-off errors have gradually shifted the fields away from the 
constraint (j).(p = 1. So, like in the planar case to correct for this we have 

rescaled (f ^ (f/ \J (f.cf every few iterations. Each time, just before the rescaling 
operation, we have evaluated the quantity = — 1 at each lattice point. 
Treating the maximum of the absolute value of /U as a measure of the numerical 
errors, we have found that max|/i| ~ 10~^. This magnitude is useful as a guide to 
determine how reliable a given numerical outcome is. Usage of unsound numerics 
in the Runge-Kutta evolution shows itself as a rapid growth of max|/i|; this also 
occurs, for instance, when the solitons pinch-off. 

The parameter A in © has been set to A = (1, 0) all throughout. 



4 Three solitons 

First we have considered states with three solitons. Our initial configuration is 
given by taking k = 3 in the elliptic function of order 3 

^^ aiz-a,)aiz-a,)a{z-as) ^ (12) 
aiz - 61) aiz - 62) c7{z - 63) f^^ J^i 

The values of aj , hj have been selected in such a way that the solitons lie sym- 
metrically around a circle in the network (|11|) . This is easily achieved by fixing 
ai,6i to reasonable values and setting 



(13) 



ai = ai — c \ 
b[ = 61 - c J ' 

where c = (2, 2) stands for the centre of the period cell. Then 

This symmetrical arrangement gives solitons of the same size, and satisfies the 
selection rule in 1)121) for any choice of the complex numbers a'l, 6^ 

3 



= j:b'r (15) 



K = l 

3 



ft=l 
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Next we go back to a,- and bj through 



= + C ^gx 
bj = b'j+C ^^^> 



and supply the system with an initial speed vq by boosting 

Vj = -voey:p{i(3j). (17) 



Qj — > aj + Vjt 
bj bj+Vjt 



It is now possible to evaluate the time derivative of W at the initial time. In- 
serting ((T7|l into (fT^ we get 

3.W(t) Uo= - n .,1^^;^^ - "-'^^^W,. (18) 

(T{z-aj) a{z-bj) 

where Wj = a{z — aj)/a{z — bj) denotes the j'-th soliton. Our initial-value 
problem is defined by (fT^ and ((TH|l . 

Choosing 

ai = (3,2) ] 

6i = (1.45,1.95) I, (19) 
vq = 0.35 J 

gives an initial configuration whose energy density is exhibited in figure ^ The 
solitons are placed on the vertices of an equilateral triangle; the first lump is 
situated along the line y=2 at /3i =0°, with respect to which solitons 2 and 3 are 
rotated j32 = 120° and /^a = 240°, respectively. These angles are readily checked 
from the picture with the help of a protractor. 

The results of simulations for this case are depicted in figure HI for various 
times. The three solitons are fired towards the centre of the triangle and in so 
doing they expand: the peak Emax of the energy density decreases. The solitonic 
trio collides head-on and coalesces in a ringish structure, then emerging towards 
the vertices of the dual triangle, that is, the initial line of approach of a given 
incoming soliton forms an angle of 7r/3 with the line along which an outgoing, 
emerging soliton progresses. This 7r/3 scattering can be best appreciated in figure 
121 where both the initial state (t=0) and the final state {t = 2) are displayed 
together. 

Returning to figure |21 we see that Emax becomes narrower with time, par- 
ticularly after the lumps scatter off and start drawing away from each other. 
At t = 2, for instance, the maximum value of the energy density goes up to 
Emax = 579.37. Soon after t > 2 this peak gets so spiky that the numerics 
breaks down: the instability of the planar 0(3) model takes over and leads to 
singularity formation. 

Our result is noteworthy: the initial configuration, although positioned at 
the vertices of an equilateral triangle in the period cell, does not produce D^, 
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symmetry because the torus itself, being homogeneous but not isotropic, has no 
such symmetry (the fundamental grid has directed sides). One could reasonably 
expect the lumps to scatter along directions that need not respect D3 symmetry. 
Therefore, the rationale applied to explain tt/N scattering for the model on ^2 
is no longer valid for the model with periodic boundary conditions. 

However, we observe that as the solitons are well localised the boundary 
conditions may not be very important. Note that our numerical results are 
consistent with this expectation. To test this further we could place our solitons 
nearer the edges of the grid and see whether we still observe the 60'^ scattering. 
We hope to investigate this issue in the future. 



5 Four solitons 

Next we have looked at the = 4 configurations. The initial field is given by 
the elliptic function of order 4 

a{z - ai) a{z - 02) ajz - 03) a{z - 04) A A 

where aj , bj are chosen so that the solitons sit symmetrically at the vertices of 
a square in the basic cell. A treatment parallel to that of the previous section 
leads to 

aj = «lexp(iv.,)l = = 1,2,3, 4, (21) 

The condition between the zeros and poles is again verified V a'l, b'l : 

Ea;- = (0,0) 

4 

= Eb'r 

i=i 

Both the initial velocity and the time dependence are introduced by boosting: 
^' Z l]Xlft }' ^i = -^oexp(i^,), (22) 



where equation (|16|) should be kept in mind. 

As pointed out in section |2J four lumps involve a large energy and while 
evolving in time they become too spiky before we can learn anything much about 
the scattering process. We have therefore studied this system in the Skyrme 
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version of the theory, where the sohtons are stable and may be examined for as 
long as required. 

Instead of the Lagrangian density (jlj) we have thus taken: 



C = il-U Ll-lJ L-lJ sOa- — ^ — —( WA^ + W-A^ - VF? 



(23) 



The configuration (j^J is no longer an exact solution of the field equation derived 
from ()23() . albeit it is a very good approximation to it. We should also stress 
that the presence of a Skyrme term does not affect the trajectory of the lumps 
before the collision or their scattering angle. 



(24) 



Taking 

ai = (2.05,2.55) 
bi = (2.95,2.95) 
01 = 10-3 
vq = 0.35 exp(i7r/4) 

entails the state of identical 'baby skyrmions' illustrated in figure H) The solitons 
W2, W3 and W4 are rotated 90°, 180° and 270° with respect to Wi, which we 
have conveniently placed in the first quadrant, roughly on the central diagonal 
joining the grid points (0,0)-(4,4). 

The system (|2()|) gets moving via (|22|) and zeroes in on the middle of the 
mesh, where the four skyrmions bump head-on into each other. This dynamics 
makes the solitons scatter off and emerge towards the vertices of the dual square, 
emerging at 45° with respect to the initial direction of motion, as depicted in 
diagram 13 Figure El superimposes both the initial state (t=0) and the final state 
(t = 3), allowing greater clarity in the appreciation of the scattering angle 7r/4. 

Unlike the outcome of the previous section, the dual square scattering on T2 
is not so surprising. For although the initial data doesn't have D4 symmetry, it 
does have 4-fold rotational symmetry and zero angular momentum. Note that 
the boundary conditions Q break the 5*0(2) rotational symmetry of the plane 
into a 4-fold rotational symmetry. 



6 Geodesic approximation 

Note that, in analogy to the S2 case, when v in (17) and (22) is set to zero 
our initial configuration is a static solution of the equations of motion. If we 
change aj, bj to a new value given by (17) or (22) for a particular value of t the 
new configuration is again a static solution of the equations of motion. Hence 
as V changes the changes (17) and (22) connect configurations which correspond 
to static solutions of the equations of motion. Thus it is reasonable to expect 
that a system set off with a small v will follow such a change. This expectation 
goes under the name of geodesic approximation (the system evolves by changing 
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its zero- mode parameters). Its validity is not expected to depend too much on 
whether the model is defined on 5*2 or T2. 

So far we have been concerned with solitons of equal size. Let us now look 
into the more general and interesting situation of energy lumps of different sizes, 
illustrating the proceedings by studying the case = 3. 

First we set up the initial configuration and evolve it through our standard 
numerical simulation. Then we choose a set of collective coordinates to reproduce 
the results of our simulations (trajectory, scattering), offering an explanation in 
the framework of the geodesic approximation. 



6.1 Numerical simulation 

Our 3-soliton system is still given by a function of the form H12() . but with a 
layout not as symmetrical as before. Instead of (|14|) we put 

t ; ",Mi-i)f (25) 

with j=l,2,3 and the numbers a'l, b'l being fixed as customary. The initial three 
soliton configuration can be written as 



W 



a{z — oi) a{z — 02) cr{z — 03) 
a{z - b) a{z - 62) o-{z - 63) 
a{z' - a[) a{z' - a[e'''^) a{z' - a'^e''''^) 
a{z' - b') a{z' - 6;e*"2) a{z' - 6;e-*°2) ' 



(26) 



where 



27T 

b' = a'l + 2{a[ - b[) cos ^2 [a2 = — - ^] (27) 

ensures that the zeros and poles in (|26() comply with the constraint in ®. As 
usual we switch between primed and unprimed numbers through formula (|16)) . 
plus b' = b — c and z' = z — c. We have also used the fact that exp(za3) = 
exp(— ia2)- 

Note that equating ^ to zero simplifies aj to {j — l)27r/3 = f3j = of (|14|). 
whereupon the elliptic function reduces to the field l(T^ simply because 
cosa2 = COS/32 = — 1/2 — > 6' = according to 1)27(1 . 

An angle ^ 7^ generates solitons of different sizes which no longer enjoy the 
positional symmetry boasted by the arrangement ()14j) . In the set-up ()14|) . where 
^=0 and the solitons have the same size, the required relationship X^^j = ^^'j 
looks after itself as shown in (fT5|) . For ^ 7^ a demand of the kind (|27j) is needed. 

Now we have evolved the solitons H26|) in the pure version of the CP^ model, 
with the solitons being sent into collision in the regular fashion: 

aj —>■ {a[ — vot) exp {iaj) + c 

b ^ {b'-vot) + c }; i = l,2,3; s = 2,3. (28) 

bs {b'l - Vot) exp (ias) + c 
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The associated energy distribution is shown for a typical case by the contour 
plot of figure [7| which corresponds to a choice of parameters (|19() in addition to 
C = 10°. 

The starting configuration is represented by the three wider structures whereas 
the final, scattered state is the narrower trio plotted therein. The sense of mo- 
tion is clearly indicated by the arrows. We note that for both configurations the 
solitons are rotated an extra angle ^ with respect to Wi, as compared to their 
partners of figure|3J And, unlike the latter, it is not the case that the three initial 
lumps have the same size; nor they are situated at the vertices of an equilat- 
eral triangle. Apparent as well from our simulations is that the scattering angle 
differs from 7r/3. This is a consequence of considering nonsymmetrical solitons, 
whose collisions are not elastic and thus involve energy transfer. 

Clearly the distance between the solitons differ for different pairs of them 
and so their interactions are not the same resulting in a different scattering 
pattern. Can we explain this difference? In the next subsection we show that 
an explanation can be provided in terms of, appropriately chosen, collective 
coordinates and the motion following appropriate geodesies. 



6.2 Collective coordinates 

In order to make a wise selection of collective coordinates let us consider closely 
the positions and sizes of the solitons defined in the previos subsection. 
The location of the solitons (|26() are clearly given by 

a[ + b' a[ + b[j^^ a[ + b[^_,^^^ ^^^^ 



and their sizes are 



Put 



2 ' 2 '2 
a'i-6' ,a[-b[ .a[-b[ 



(30) 



a'i = k-x, b[ = k + x (31) 
so the positions (|29jl adopt the form 

fc-x[l + 2cos(a2)], A:e^"2, ke-''''^ , (32) 

while the sizes (|5n)) read 

|2xcos(q2)|, IxI, IxI- (33) 

Consistently, for ^ = the description corresponds to objects of the 

same size \x\ which are symmetrically situated at k and A; exp(ibi27r/3) (with 
respect to the centre). For simplicity, we have taken two lumps with the same 
size and a third lump with a different size. 

A possible collective coordinate description involves treating k, x ^-nd 02 as 
collective coordinates. Thus, in the simulation we expect 02 and x to remain 
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approximately constant and k to vary. We are suggesting that the scattering 
can be understood as proceeding (on average) with only k depending on time 
and varying from (take k real for simplicity) A; > to A; < 0. It is easy to 
show, although a bit tedious in practice (this involves estimating various elliptic 
integrals which can be done partially analytically, partially numerically), that 
both on ^2 or on our torus the kinetic energy of the motion just described is 
finite. So such a behaviour is possible. 

We have carried out our collective coordinates "motion" using the specific 
values 

X = (0.55,0), fc = (n,0), [1,-1], (34) 

where u varies across the interval in steps of 0.2. We see from (|32j) and (|34|) that 
the configurations for u > {u < 0) correspond to incoming (outgoing/scattered) 
lumps. The value u ~ represents the situation where the solitons have collided 
and are on top of each other, coalescing in the centre of the mesh. It is largely 
the behaviour of at n ~ what determines the scattering angle. 

The state for A;=(1,0) (incoming lumps) and A;=(-0.2,0) ("scattered" lumps) 
are depicted in figure |H1 for the case ^ = 10°. As long as the path followed by the 
lumps and the scattering angle is concerned, the similarity with the numerically 
evolved situation of figure[7|is clear. We stress that since we are mostly interested 
in the relative position of the solitons and their scattering angles, it is immaterial 
that the breadths of the solitons in figure [7| differ from the solitons of figure |H1 

A more detailed comparison can be made by viewing a trajectory plot, one 
showing consecutive snapshots [corresponding to /c=(l,0), (0. 8,0),. ..,(-0,2), (-0.4,0). ..(- 
1,0)] of lump positions as they approach each other and scatter off. 

This is shown in figure IHl where the full collective coordinate motion has 
been ticked with small circles and the motion according to the time evolution of 
subsection 6.1 has been sketched with solid lines. Note that after scattering the 
path of the numerically evolved lumps (continuos lines) cannot be followed much 
farther; this is because the solitons get very spiky and the simulation breaks 
down. 

Our simulations have been carried out for several values of ^ and have been 
compared with the corresponding collective coordinate motion. We have always 
found the trajectories from both approaches to be in agreement, thus supporting 
our choice of collective coordinates. Another illustration is provided by figure EH 
for the case ^ = 20°. 
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Emax=97.33; t=0 




Figure 1: The energy distribution for A'^=3 at the initial time, both in three dimensional 
and contour plot forms. 
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t=0.5 



t=1 




t=1.5 



t=2 









1 







# 





Figure 2: The evolution of the three lumps of figure Q at various times. 
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Scattering at 60 degrees 




Figure 3: Scattering at 7r/3 radians. 
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Emax=123.014;t=0 



Figure 4: The energy distribution for N—A ai t — 0. 
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Figure 5: The evolution of the four lumps of figure IH gets underway. 
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Scattering at 45 degrees 



2 




2 



Figure 6: Scattering at 7r/4 radians. 
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Scattering of non-symmetrical lumps 



2 




2 

Att=0: a1=3+2i; b=1 .9397+1 .96581 



Figure 7: Numerical simulation for solitons rotated an angle ^ = 10° with respect to 
the symmetrical configuration of figure ^ They do not scatter at 60°. 
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Collective coordinates scattering 




2 



Figure 8: "Scattering" according to our collective coordinates view for ^ = 10°, to 
be compared with the numerical result shown in figure [7| The incoming lumps are 
pictured for k = (1,0) and the outgoing structures correspond to k = (—0.2,0). It is 
not important that the widths of these lumps are not the same as their sibblings of 
figure [7| we worry about the relative positions and scattering angle. 
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Trajectories 
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Figure 9: Trajectory plots = 10°) sketching the path followed by the humps pre- 
sented in figures [7| and |H1 The itinerary according to our geodesic approximation (cir- 
cles) shows very good agreement with the motion obtained via numerical simulation 
(sohd lines). 
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Figure 10: Plots comparing the trajectory obtained via numerical simulations (solid 
lines) and the trajectory gotten using the geodesic approximation (circles) for ^ = 20°. 
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7 Conclusions 



In this article we have studied head-on cohisions between soUtons in the (2-|-l)-D 
CP^ model with periodic boundary conditions, that is, with the model defined 
on a flat torus T2. 

Through numerical simulations we have found that solitons of equal sizes (at 
the initial time) scatter at an angle ir/N, or dual-polygon scattering, where 
is the soliton number or topological charge of the system. In this paper we have 
focused on = 3, 4 (the case N = 2 has been considered previously ^ and it 
has been found to comply with the 7r/2 scattering). 

Unlike the usual model on 5R2, our model with periodic boundary conditions 
breaks the SO (2) rotational invariance of the plane, leaving us with a numerical 
mesh with directed sides where the initial soliton configuration has no symmetry 
under the dihedral group -Dat. So it is remarkable to still find dual-triangle scat- 
tering in this format. On the other hand, since the 5*0(2) symmetry breaks into 
a four-fold rotational symmetry, the dual-square scattering is less unexpected. 

We have also considered solitons of different sizes (at the initial time) and 
have observed that the scattering angle in no longer n/N, outcome arising from 
the fact that there is energy transfer in collisions between unsymmetrical solitons. 

By reparametrising the quantities describing the positions of the solitons using 
a juditious set of collective coordinates, we have been able to reproduce the above 
numerical results, thus offering an explanation of the scattering process. We have 
illustrated this approach using a 3-soliton field. 

These results raise important questions pertaining to the the interplay be- 
tween the symmetry of the initial configuration and the lack of symmetry of 
the torus itself. As pointed out at the end of section 0J the non-isotropy of the 
torus might affect the evolution of the lumps if they are initially placed near the 
boundary of the mesh. What about systems with A^=5,6,...? Numerical experi- 
ments on a periodic, rectangular grid would also be worth performing. We hope 
to report on these matters in the near future. 
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